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gl.  Introduction. 

A  large  structure  is  formed  by  the  coupling  of  simple  structural  elements.  In 
thia  paper  we  study  the  simplest  type  of  such  structures  which  is  made  up  to  two 
coupled  strings  modelled  by  quasilinear  or  linear  wave  equations.  We  install  two 
stabilizers:  one  at  the  (left)  boundary  and  one  at  an  in-span  point.  We  wish  to  study 
the  exponential  stability  property  of  this  coupled  dynamic  structure.  The  nonlinear 
or  linear  partial  differential  equationa  are  described  below: 


(nonlinear) 


(«.t) 


a* 

liT 


(x.t) 


0  X  ^  1, 


1  <  X  ^  2, 


(1.1) 


(linear) 


-  cj  =  0.  0 .  X  ^  1. 

f^<x.t)  -  cj  =  0,  1 X  X  ^  2. 


(1.2) 


where 

o,,  Oj  satisfy  ei(0)  =  0,  el(u)  ^  0,  i  =  1,2, 

and 

c,  =  VTTTp  ,  c,  =  vT^Tp 

T|  =  tension  constant  on  string  i,  i  =  1,2, 
fi  =  maaa  density  per  unit  length. 


Supported  in  part  by  NSF  Grant  84-01297A01  and  APOSR  Grant  8S-02S3.  The  U.S. 
Government  is  authorized  to  reproduce  and  distribute  reprints  for  governmental 
purposes  notwithstanding  any  copyright  notation  thereon. 
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Wo  hove  UiciLly  uuaumod  Ihol  Iho  muRO  donaitiou  on  both  uLringH  ore  iduriLicul.  ('t'hiR 
information  is  contained  in  oi  and  o,  in  (1.1)).  The  length  of  each  string  has 
been  normalized  to  1. 

At  the  left  end  x  =  0,  a  stabilizer  is  installed  satisfying,  respectively,  the 
following  condition 


<t.(w,(0,l))  -  l<>t(0,l)  =  0  (for  (1.1)),  kj  I.  0, 
cJwm(0,1)  -  kJwjCO.t)  =  0  (for  (1.2)),  kj  >•  0. 


At  the  intermediate  node  x  =  I,  another  stabilizer  is  installed  according  to  one  of 
the  following  two  seta  of  dissipative  transmission  conditions  (compare  (2))  hold: 


ffi(w.(l-,t))  =  e,(w.(l*,t)). 

wi(l-,t)-wt(l\t)  =  -kJo,(w.(l-,t))(=  -kJo,(w.(l+,t))),  kj  *  0, 


w(l-,l))  =  w(lM), 

Oi(w.(l-,t))-o,(w,(l+,l))  =  -kjwt(l-,l)(=  -k!wt(l+,t)),  kj  »  0, 


for  the  nonlinear  system  (1.1).  For  the  linear  system,  the  counterparts  are 


clw,(l",t)  =  cjw,(l'*’,t) 

wi(l-,t)-wt(l  +  ,t)  4  -kjc?w.(l-,t)(=  -kicjw.d+.t)),  kf  s  0, 


w(l-,t)  s  w(l*,t) 

c!w.(l-.t)  -  c;w.(l\t)  =  -klwt(r,t)(=  -kjwtd^,t)),  kj  s  0. 


Note  that  if  k|  =  0  in  (1.5),  (1.5)’,  (1.6)  or  (1.6)’,  then  there  is  no  loss  of 
energy  at  x  =  1  and  the  joint  is  conservative.  We  call  ko  and  k'  the  feedback 
gains  at  x  z  0  and  x  z  I. 

At  the  right  end  x  z  2,  we  may  assume  that  it  is  either  fixed  or  free: 


w(2,t)  =  0  (fixed  end) 

w,(2,t)  =  0  (free  end) 


We  want  to  study  the  effects  of  stabilizers  for  coupled  nonlinear  and  linear 
vibrating  strings  as  described  above.  We  will  be  primarily  concerned  with  the 
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fullowing  question: 

IQJ  "When  does  the  solution  of  the  shove  coupled  wove  equations  decay 
exponentially?" 

The  answer  tx>  the  above  question  is  well  understood  in  the  cae  of  a  single 
(nonlinear  or  linear)  vibrating  string.  Greenberg  and  Li  [6]  proved  that  for  (l.l), 
(1.3)  and  a  conservative  boundary  condition  at  x  =  1,  the  solution  decays 
exponentially  in  the  C'-nora  for  any  k'  ^  0,  provided  that  the  initial  data  is 
sufficiently  small,  smooth,  and  that  appropriate  compatibility  conditions  are  satisfied. 
In  the  linear  case,  uniform  exponential  decay  of  energy  follows  as  a  simple  exercise 
of  the  method  of  characteristics. 

The  problem  of  coupled  vibrating  strings  is  equivalent  to  that  of  a  hyperbolic 
system.  Some  recent  papers  by  Chen-Coleman-West  [2]  and  Qin  [91  have  provided 
partial  answers  to  [Q]  under  study  here.  Another  recent  paper  [7]  by  K.L.  Huang 
has  also  introduced  a  direct  way  of  proving  exponential  decay  of  solutions  for  linear 
systems.  His  method  is  to  establish  a  uniform  bound  of  the  resolvent  operator  on  the 
imaginary  axis.  In  this  paper,  we  will  provide  answers  to  [Q]  along  the  directions  of 
[2],  [7]  and  [9]. 

The  organization  of  our  paper  is  as  follows.  In  §11,  we  first  give  a 
counterexample  of  Qin's  theorem  in  [9]  and  state  Its  correct  version.  In  illl,  we 
apply  the  corrected  theorem  of  Qin  to  coupled  nonlinear  vibrating  strings  and  state 
soma  sufficient  conditions  for  exponential  decay  of  solutions  with  small  data.  In  §1V, 
we  apply  Huang's  theorem  to  givo  a  complete  answer  to  [Q]  for  coupled  linear 
strings. 


§11.  A  Counterexample  and  Correction  to  a  Theorem  by  T.H.  Qin  [9]  for  Qiiasilincnr 
Hyperbolic  Systems  with  Disainative  Boundary  Conditions 
Let 


+  A(u(x,t))*“^*‘‘'^^  =0,  0  X  X  1,  I  0  (2.1) 

O  t  VK 

be  a  final  ord«r  quoailiiionr  hyperbolic  Ryatom,  whoio  ii  c  9i**  A{«i|  in  iin  N«ri 

aufficienlly  ainooth  niulrix  forictiori  of  the  vuruibin  u  only. 

As  in  (9K  aaaunie  that 

(Al)  Syalem  (2.1)  ia  hyperbolic  for  aufficienlly  ainiill  lul  in  (he  fnilowintf  Horim': 
i)  The  matrix  A(ti)  hee  N  RnuMith  roal  oi^onvalucH  ^  nrxl 


A,(0)  4  ...  4  A«((n  /  0  ^  A^»,(0)  4  ...  4  A^(0). 


/  X .  2  1 


Nn/v'i ■V'". ^  ^ ^ 


•H#  -VvA-  .• 


:‘‘\,k  k-  .  -i  -r-  •■: ' 


ii)  A(u)  hu8  N  lirioariy  indopendant  iefl  eiKonvectors 
^j(u)  =  ((ji(u),...,fj„(u)),  1  *  j  4  N 

corresponding  to  each  real  eigenvalue  Aj(u). 

Without  loss  of  generality,  we  asBuma  that  the  matrix 


<,.(U). 

. f.H(u) 

. (mm(u) 

u  =  0, 

i.o., 

'  1 

“1 

U| 

u»+i 

• 

.  u"  a 

Uai 

ua 

Dua  U>  (2.2)  and  (2.5),  ganaral  boundary  conditlona  for  syslem  (2.1)  should  have  the 
following  forma: 


u*'  s  F(u')  at  K  -  0, 
u‘  =  C(u«<)  at  X  =  1, 


K  :  R"  ^  R"- 


G  :  R"  "  .*  R" 


('•1 

»a 

.  «  = 

tie 

are  (;‘-amof)t.h  veclor-valuod  functions.  Wo  dofino 


I’l  •  a)«(N  a)  autrix. 


<ir<5y  Ti* 'if 
C.  -'.r  \;>K-  « 
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til)  matrix, 


where 


ii_  ,  ■*(►-> . f.) 

iu*  i(ut,...,u.)  ’ 


au”  a(u.+, . um) 


Let  the  initial  condition  of  (2.1)  be 


u(x,0)  =  *<x),  0  a  X  a  1.  (2.8) 

In  addition  to  (Al),  we  aeaume 

(A2)  The  initial  and  boundary  conditions  (2.8)  and  (2.7)  satisfy  the  following: 

i)  F(0)  =  0,  G(0)  =  0. 

ii)  At  X  =  0  and  x  =  1,  the  following  compatibility  conditions  are  satisfied: 


♦  ”(0)  =  K{#*(0)).  ♦«(!)  =  G(*'>(1))  (2.9) 

[A‘*(#(0))  -  |^{♦*(0))A“(♦(0))}^1|^  + 

(2.10) 

[A”(#(0))  -  |^(»»(0))A’»(»(0))]^^°^  X  0. 


(A‘'(*(l))  -  |^(♦Il(l))A»•(♦{l))]i|^  + 

[A'»(f(0))  -  |^(#tl(l))A”(a(l)))^*^^  =  0, 


(2.11) 


where  A’’,  A”,  A”  and  A’*  are  blocks  of  sub-matrices  of  A  corresponding  to 
the  decomposition  u  =  (u'.utt): 


A(u) 


A"(u) 

A>'(u) 


A”(u) 

A”(u) 


□ 


For  any  nun  square  aatrix  M  =  (B|j),  we  define  its  absolute  value  aatrix  H 


by 


•  n^r\Mrym  W  WITU  *r\(  m#  JW  W¥irCfa: 


ft  •  ('■♦jOn«f»- 


This  mairix  M  satisfies  the  properly  that 


I.  s  I  Max  £  for  all  v  c  ft", 

j  =  i  •*  ^ 


I 

m 


Ivl,  =  Max  iv,l  for  v  =(v,,...,v„)  e  R", 

I<l4n 


We  thua  deftna  a  norm  I  - 1  for  M  by 


■Ml  ■  aax  £  (■  IMI). 

•  4(4n  j  =  «  •* 

Wa  now  stata  tha  “thaoram”  of  T.H.  Qin  in  (91. 


Theoram  0.  Aaauma  (AI)  and  (A2).  Aaauma  further  that  the  apectral  radii  of  B| 
and  Bt  are  leaa  than  1.  Then  there  axiata  <  ^  0  auch  that  tha  mixed 

initial-boundary  value  problem  (2.1),  (2.7),  (2.8)  admita  a  unique  global  amooth 
aolution  u(t,x)  for  t  *  0  and  lu( t ,  • )  1^,,  decaya  exponentially 


?l.>l 


■  u(t,  )l^,,^Q  a  Ke-«t.  K,a  X  0.  for  all  t  k  0, 


(2.12) 


provided  that  i)  *  *^'*C*(0  1}  ^  ^  depends  on  4  and 

the  rate  of  decay  a  is  at  least 


In  » 


(2.13) 


X.,,,  ■  Bin  lAtl-o/  (for  fl’  V  0  suae  Haul  1  nuvlier)  (2.11) 

ttf  I^N  ^  *■"*  •  ®  ^ 

p  is  the  smallest  positive  Integer  making 


Bax(iBfi,iBSi)  ■  a.  4  1 


and  0  is  any  real  number  in  {ootl)» 


11  is  well-known  thot  solutions  to  quasilinear  hyperbolic  systems  with  smooth 


I 
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data  develop  shocka  within  finite  duration.  Therefore  it  ie  required  in  the  cunditiona 
ot  Qin's  theorem  that  the  C‘-norm  of  #  be  sufficiently  small.  I'hus  normally 
exponential  stability  (2.12)  can  only  be  expected  for  small  data. 

The  condiliona  on  the  spectral  radii  of  B|  and  Bj  in  the  statement  of 
"Theorem  0"  seem  so  intuitive  and  natural,  as  we  anticipate  that  waves  would  lose 
strength  exponentially  after  repeated  reflections  on  energy  absorbing  boundaries. 
Unfortunately,  the  conclusion  is  erroneous  due  to  the  fact  that  in  a  system  there  are 
several  waves  travelling  with  different  speeds,  thus  their  superpositions  may  form  a 
wave  which  is  undamped.  The  following  is  a  counterexample. 

We  consider  a  linear  system  formed  by  the  coupled  wave  equations  (cf.  (1.2)) 


i*w(x.t)  ad^w(x.t) 

)t>  ■  ax>  -  “■ 

a^w(x. t)  _ 

it'  Jk' 

(cj  c\) 


0  ^  X  ^  1 ,  t  0 
I  ^  X  2,  i  X  0 


with  boundary  and  tranamiaalon  condiliona 


(2.16) 


w(0,l)  =  0 
w„(2.t)  =  0 

=  -kjcjw.d-.t) 
cfw,(l“,t)  =  cjw,(l^,t) 


at  left  end  x  «  0, 
at  right  end  x  =  2, 

at  X  =  1,  cf.  (1.6). 


(2.17) 


Let  us  IrunsforDi  them  into  a  hyperbolic  aystem  by  letting,  for  x  c  (0,1), 


«i(K.t)  =  i|-cj|^(2-x,t)  ^  |Y>i(2-x,t)j  , 

“i(x.t)  *  ft*'*’*’ '■^1  > 

w>(x.l)  =  *  ft*'**’ ■ 

W4(x,t)  =  |[‘=»f^*2''‘.f)  +  f^(2-x,t)|. 

Then  v»  :  (w,,W],w,,W4)  :  (w*,*”)  (w*  =  (w,,w,),  w>l  =  (w,,w.))  aaliafieo 


v(x,t)  + 


f-c,  0| 


C|  Jx 


((X,t)  =  0 


w*'(0.t)  =  [°  *]w«(0,t)  •  Dow*(0,t) 


=  (c,+c,+kjc,ca)  *1 

■  D,w"(l,t). 


2ct  -c,+Ci+k,C|C; 

c,-cj+kfc,cj  2cj 


It  is  easy  to  see  that  (Al)  and  (A2)ii)  are  aatia/ied.  (A2)(it)  is  satisfied  if  the  initial 
condition  is  chosen  to  satisfy  the  compatibility  conditions.  (2.5)  is  trivially  satisfied. 
NoW|  DuO|  is  cooiputcd  to  be 


=  DoD,  =  -(c,+c,+k,c,c,)  ‘ 


c.-Cj-tkiCiC,  +2cj 


c,-c,-k,c,c. 


which  haa  aigenvaluaa 


=  -(ci«-c,*kJc,cj)'*<c,-c,»[c,c,(ktciC,-4)]‘/').  (2.20) 

Bj  :  0|Da  haa  identical  ei^anvaluea  aa  B|  :  DaD|  aa  given  above.  It  la  eaay  to 
aee  that 


*  1.  for  any  k,  .  0. 


Hence  all  of  the  hypothaaea  of  "Theorem  0"  are  aatiafied. 

We  now  prove  that  aolutiona  to  (2.18)  do  not  decay  exponentiully  by  ahowing 
that  the  coupled  wave  syatem  (2.15),  (2.17)  haa  at  loaat  Home  eigonvaluea  located  on 
the  ioiugiiiary  axie.  Wo  write  an  eigenaolulion  of  (2.16),  (2.17)  aa 


w(x.t)  = 


e^‘#,(K),  0  a  X  «  1 

e^*#a(x),  1  a  X  *  2. 


♦,(x)  -  (t/c,)«4,(x)  =0,  0  .  X  a  1, 

♦»(x)  -  (A/c,)'4,(x)  =0,  1  a  X  a  2. 


Due  to  the  firat  two  boundary  conditiona  in  (2.17),  we  have 


Wi:vw>!»y.H«T>i 


S’ 


1<) 


If 


£i  -  i  -  2n 
Cj  "  M  2n+l’ 


(L  =  2n,  M  =  2n+l) 


then  M+L  =  4n+l,  4-L 


1.  c,+Ca 


4n>  1 


2n<-l‘^'’ 


and  (2.28)  gives 


(4n+l)coa(4n+l)Q  -  cosO  =  0 
ain(4n'>l)D  -  sinO  =  0. 


(2.29) 


It  is  obvious  that  0  =  ”/2  is  a  solution  of  (2.29).  Hence 


A  =  =  i(n  +  1)"^, ,  any  positive  integer  n, 

is  an  eigenfrequency  of  (2.16)t  (2.17)  which  is  undamped  as  He  A  =  0.  This  is  a 
counterexample  to  "Theorem  0”. 

The  fallacy  of  "Theorem  0"  is  due  to  the  fact  that  in  [9]  the  author  has 
forgotten  to  take  certain  absolute  values  between  steps  (52)  -  (54)  in  [9i  pp.  295- 
296].  The  aistakes  can  be  easily  corrected  by  changing  Bi,  Ba  in  the  atateaent  of 
“Theorea  0"  to  fii,  Ba.  He  state  the  corrected  version  of  Qin’s  theorea  below: 


I 

i 


Theorea  1 .  Asauae  (Al)  and  (A2).  Assuae  further  that  the  spectral  radii  of  B|  and 

S,  are  leas  than  1.  Then  there  exists  4^0  such  that  the  aixed  initial-boundary 

value  prublea  (2.1),  (2.7),  (2.8)  adaits  a  unique  global  aaooth  solution  u(t,x)  for 

t  k  0  and  lu(t,')l-i  decays  exponentially 

b  (0,1) 


lu(t,  )ij,,^g  4  Ke-“».  K,a  k  0  for  all  tkO, 


(2.30) 


provided  that  •♦•(...(o  i)  ^  '♦’'c®(0  1)  ^  *’  (2.30).  k  depen.ls  on  4  and 

the  rate  of  decay  a  is  at  least 


In  <r 


2p 


with 


■  ,5j2„'*''c®(-9’,a’)  ’ 


(for  Borne  saoll  4'  ^  0) 


p  is  the  sauV.lest  positive  integer  nuking 


(2.31) 


1 


t 


aax(IBfl'li?l)  -  a.  1 


(2.32) 


00 


and  a  is  any  real  number  iti  (OoiD*  □ 

Note  that  the  above  correction  imposes  very  severe  restrictions  on  matrices  B| 
and  Bj.  It  does  not  apply  to  many  cases  studied  in  [2]  e.g.  as  one  realises  that  in 
[2]  the  spectral  radii  of  B|  and  Bj  arc  less  than  1,  but  those  of  B|  or  §,  are 
usually  not  less  than  1. 


9111.  Coupled  Nonlinear  Vibrating  Strings  with  Point  Stabilizers. 

Let  us  first  consider  the  system  (l.l)i  (K3),  (1.5)  and  (1.7).  We  will  transform 
the  system  into  a  form  which  allows  the  application  of  Theorem  1. 

Without  loss  of  generality^  we  assume  that 

a|(0)  k  oi(0).  (3.1) 


y,(x.t)  =  -  ^(2-x.t). 

(3.2) 

z,(x,l)  =  z,(x.t)  =  |^(2-x,t), 

for  X  e  (0,1),  t  ^  0.  Uaa  tho  following  Riemann  invariants  to  diagonaliza  tha  aystam: 


Ut  •  5(*i  +  J  '■'oTTSTd’j). 

0 

y< 

u»-)  •  |(Z.  -  I  '/STT^dn), 


i  =  1,2. 


Then  u  -  (u,  ,U]  ,U),u«)  :  (u*,u>*){u’  :  (u,,Uj),  u*'  :  (u,,u,))  natlufie 


-c,(u,-u.)  0 

0 

0 

0  -c,(u,-u>) 

0  0 

0 

0 

^  n 

Ci(u,-u,) 

r — u  =  0 , 
ix 

0  0 

0  ^  X  ><  1,  t  i  0. 

0 

c, (u, -u, ) 

c,  (y)  =  v'e?yi 


i  =  1.2, 


r.'.v.vjv.v 


.  .iji^.’!.:  ^  A'  :^r 


*  ■'  '  '  '‘K*!  .'  '•  ;,'  ■  ■, 

'4 ''■'-^:..,v,.  -  ■  ■  ' 


v'-  ■■•— 


✓ 

'jf 

'j 

» 


denotes  the  wave  speed  on  slrinfi  i,  and  U(v)  is  defindod  implicitly  by 


u(y) 

V  =  J  Vo|  (7))  dj). 


We  now  transform  boundary  and  transmission  conditions  (1.3),  (1.7)  and  (1.5)  in 
terms  of  u: 


Uj(0. t) 

u"(0.t)  =:  =  F(u*(0,t)) 

u*(0.t)J 


fU|(*.t) 

u‘(0,t)  =  =  G(u”(0,t))  = 

lu,(l.l)J 


F. (u‘(0,t)  1 

K.(u»(0.t))J 

G, (u”(l.t) 
0.(u"(l.t)) 


where  F(0)  =  G(0)  r  0,  and 


F,(ut)  X  F,(u,)  X  -u. 


iu.  “*  Ju, 
F4(u«)  =  F4(u,) 


ilx  ,  £i.lUj-F,(Ut»-.»)i  (liF.!  ^  1 

jU|  c,(u,-F4{u, ))*!(*  '  ll«U|l  M’ 


JGi  _  _ 2c  a _ 

ill,  C|+Cj*l(|c,Cj 

JG,  _  C|-Ca*)l|C|Ca 
JU4  C,  ICj  *l<'c,Cj 


JGa  _  -Cl  *Ca-*)(iCiCa 
jU]  C,  tC,*ll'c,Ci 

JG  a  _  _ 2ci _ 

Hu.  ■  C,+Cj*l(*C,C, 


<-i  -  <'t  (G,  (u‘  '  )  -II,) 

Cj  =  Cj(G,{u'')-u,). 


Therefore 


B, 


[c,+Ci«k,c,Ca)  '  *  o 


c  +k* 

I  o 


Ic, -k J I 


lc,-c,-k,c,cjl 


K.K,  2cj4-‘ 

2c.4“‘K,  Kj 


(3.8) 


where  in  the  above 


A  ■  c,+Cj+k,c,Cj 

K,  »  A'‘(c,-ci+kjc,cj) 

(3.9) 

K,  •  A“* I c,-Cj-k,c,Cj i 
K,  •  lc,-kj i/(c,+kj). 

The  spectral  radii  of  Ht  and  Ba  are  determined  from  lAI  of 

det(A-i,)  =  det(A-ij)  =  A'-(K, K,+K , )A4K,KjK,-4c, c, A->K,  =  0.  (3.10) 

According  to  Theorem  1,  a  sufficient  condition  for  exponential  decay  of  solutions  with 
amall  data  ia  that  the  roots  A,,  A,  of  (3.10)  satisfy 


max( I  A, I , I  A] I )  ^  1. 


(3.11) 


The  above  can  be  determined  by  the  stondard  Routh  stability  criterion  in  automatic 
control  as  follows.  Let  us  use  the  following  Mobius  transformation  mapping  the 
interior  of  the  unit  circle  into  the  left  half  plane: 


5-1 


Substiluting  (3.12)  into  (3.10),  wo  get 
AqZ*  +  A|2  +  Aj  =  0 


(3. 13) 


where 


Ao  ■  UK|K,  +  K,+K,KaK,-4c|CjA"*K5  0) 

A,  ■  2(l-K,K,K,»4o,c,A~*K,) 

A,  ■  l-K|Ks-K,>K,K,K,-4c,c,A"*Ks. 


(3.14) 


Condition  (3.11)  ia  now  equivalent  to  that  (3.13)  has  no  rcx>t8  on  the  open  right  half 
plane.  Thie  can  be  delerniined  by  the  Routh  criterion  aa  follows: 

Ao  Ai 

A,  0  U,  ■  **  *;  *»  °  =  A,  (3.15) 

A. 

Bj _ 0_ 


Therefora  (3.11)  ia  aatiafiecl  if  and  only  if  tho  following  condition  ia  aaliafied; 

The  first  colunn  A.,  A|,  B,  {-  A,)  in  (3.15)  ia  of  one  aign 

and  A,  *  0,  A.  *  0,  i.e.,  all  of  coefficients  A.,  A|  and  A.  (3.16) 

in  (3.13)  are  nonzero  and  positive. 

We  aummarizs  the  above  in 


Theorea  2.  Assune  (Al)  and  (A2)nsin$ll.  If  l<’»0  and  k’  *  0  in  (1.3)  and 

-  O  I 

(1.5)  are  chosen  such  that  (3.16)  is  Botisfied,  then  the  oxponefilial  decay  properly 
(2.30)  h(jldH  f<ir  (t.l),  (1-3),  (1.f>)  Mml  (1.7)  provjdf^H  tluit  Mm  iriilinl  dnln  #(x) 
Hntisfioa  coinpntibilit  y  corulitionH  niwl  !«mm  Htirfifient  ly  HtiiMil  «“-nrMtn  u 

Kr  i»m  I  t  I  fill  (3.10)  it  ft  not  >1 1  t  I’  1 1  n  1  I  t  i»  t  l<  t  i  i  (  k  *  *  0 ,  w**  run 

alwnys  I  in«(  k*  '  0  surh  ([(.lO)  m  mit  i  ■;  (  i  •••I .  1  Inr  o  t  ot  •• 

exponential  stability  (for  Rfital)  data)  tori  alwiiyR  Iw'  n('hi«*veil  by  lining  a  singt** 

3 1 nh  i  1  I  7.«‘r  at  t  )»«'  I  ••  (  I  »*im|  k  O  In  i-mi  l  i  r  n  l  m  ,  w**  I’v om  i  im*  l  !•••  f  n'o •  I; '  M 


A,, 


^  0 


A,  X  I 


i'-i  '-ill  .  iiiiJilil 

«  ‘k,*  ( <'i  '  >  j  I  *  '  I  '  k'  ' 


vi:f^ .. '  ■'  -.'.•  .  -«  ■•’  • 

sjSs^^v'SW 


*  *v  i’  ' 

Eit'ij?*''  ■:  •,  :'ff«  *i' 


--Ji 


♦  UlA-lli 
ic,  ♦(:, 


iili_iiul 
c ,  ♦  k  ^ 


1*1  *<-1  kpi 

(<:,4Cj)^  c,»k' 


If  k*  is  chosen  efiuol  (or  very  close)  to  C|,  we  will  always  get  A|  ^  0,  Aj  ^  0. 
This  is  well  known  in  wave  propagation  theory  as  k*  =>  C|  in  (1.3)  corresponds  to 
the  characteristic  iapedence  boundary  condition  which  causes  Baxinuffl  energy  loss  of 
waves  at  x  =  0.  Indeed,  the  closer  k*  to  C|  is,  the  smaller  the  power  p  in 
(2.32)  becoses.  Therefore  the  rate  of  decay  a  in  (2.31)  will  become  larger.  On  the 
other  hand,  if  k^  is  not  close  to  C|,  then  A|  and  Aj  may  easily  become 
negative,  therefore  Theorem  2  will  no  longer  be  applicable.  This  is  in  sharp 
contrast  to  the  linear  case  (cf.  Theorem  7  later)  where  k*  ^  0  can  be  arbitrary 
when  k*  =0. 

I 

What  happens  if  k*  *  0?  Can  we  get  exponential  stability  by  using  only  one 
stabiliser  in  the  middle  (x  =  1)?  The  reader  can  easily  check  that  condition  (3.16) 
is  never  satisfied  no  matter  what  value  k‘  ^  0  is  chosen.  In  fact,  in  general  the 
solution  will  not  decay  exponentially  as  the  linear  counterexample  in  Sll  has 
already  shown.  One  might  wonder  whether  nonlinearity  would  work  any  differently.  The 
answer  is  still  no  as  the  following  theorem  suggests. 


Theor<im  3.  Assume  that  kj  =  0  in  (1.3).  Then  there  are  C^-solutions  to  (1.1), 
(1.3),  (1.5),  (1.7)  which  are  undamped  for  any  k*  ^  0,  for  some  nonlinear  Oi,  9%. 


Proof.  We  first  note  that  when  k*  =0.  (1.3)  beromrs 


w.(O.l)  =  0. 


We  first  (iHNiiMie  that  lh<*  f<»U<>wing  bourKhiry  vtilue  prolilc 


-r7w(x,l)  ~a|7-w(K.M|  0,  D^h''  I 

d  t  ^  «/K  I  rfx  I 

w((!,  t  )  n  t  f  N  i 

w,  {  1  .  I  )  0 

huH  fill  ('*  Hohiliori  fur  ft.  IImmi  w»*  crin  <  t  i\i\  uu.liuni 

soiuliun  for  (l.l),  (1.3)’,  (1.7)  hr  fol|f>wM.  I}‘>fit»e 


w(  K ,  I  )  . 

y( x . t  )  -  wf  X , I ) 

wf 2  V, «  )  , 


1  X  ^  0. 
0  /  X  '  1. 


I  '  K  ' 


i*i5' 

iSji* 

& 


-r-r'  ■ 

;  Vf^'vw^.;^*;'.-;  X'* -  ■ 

„  .^i  ••'  <r.‘.  *  •  . 

s-^>,%frii;'i'''';(i'''  r 


v.  -■' 

“'r  ■ 


',».*>J*’.'?WW>?I*K 


y(2x-  1 .  t ) ,  U  ^  X  ^  1 . 


1  ^  X  ^  2. 


Al  X  =  0. 


w(0-,t)  =  w(0*.t)  =  0,  w.(0-,t)  =  w.(0+,l), 

wt(0-,t)  =  vii(0*.l)  =  0,  w,.<0-,t)  =  w..(0+,t), 
wn(0-,t)  =  wi.,(0+,l)  =  0, 


BO  y  c  C*(-l,l)  and  w  c  C’(0,1).  w  aalisfies 


B  -  m  oil  1^)1  =  “■ 
0-kW-|^)]=o. 


Oxxxl,  txO 


1  X  X  2,  t  X  0. 


<»i(>l)  "  I  ff(|  »)), 

Ot(v)  •  -a(-v). 

Then  ei(0)  =0  and  a!  x  o>  i  =  1,2,  and  w  aaliafies  (].l).  At  x  :  0  and  2, 


w.(0,t)  =  2yx(-l,t)  =  2w,(l,t)  =  0. 
w(2.1)  =  y(O.t)  =  i(0,t)  =  0. 

Al  X  =  1, 

«.(!-. 1)  =  w.(l*,l)  =  0 

w,(l-.t)  =  Wt(l‘.t), 

80  (1.5)  ia  satisfied.  Hence  w  is  a  C'-soluton  of  (I.l),  (1.3)',  (1.5),  (1.7)  wliicli  is 
undamped. 

But  Greenberg  [5]  has  proved  that  with 


V  x  -A,  Co  x  0,  A  X  0, 


(3.17)  has  a  C'solulion  which  is  undamped. 
The  proof  of  Theorem  3  is  complete. 


'I  t;  - 
'  r’.‘  1  '  i.  '1  •  . 


Therefore  we  can  aafoly  staio  that  in _ genorai  nxponofiliui  ala  bi  lily _ rnnnol  be 

attained  by  uaing  only  one  atabilizer  at  the  middle  of  the  apan  for  coupled  (linear  or 
nonlinear)  wave  eguationa. 

Although  we  have  uaed  the  Riemonn  invarianta  u  in  (3.3)  to  give  proofu  end 
deduce  stability  results,  the  plain  fact  is  that  here  as  in  (6]  and  [9]  and  in  many 
other  similar  papers  the  approach  is  essentially  based  on  linearization  about  the  zero 
solution  u  :  0.  Theorem  2  is  but  one  of  the  many  examples  of  the  principle  of 
linearized  stability  which  is  valid  for  a  large  class  of  general  nonlinear  distributed 
systems. 

The  case  when  the  intermediate  conditions  and  boundary  condition  are  changed 
to  (1.5)'  and/or  (1.7)'  can  be  studied  similarly  without  any  difficulty. 


IIV.  Couplied  Linear  Vibrating  Strings  with  Point  Stabilizers. 

In  this  section  we  study  the  exponential  stability  problem  for  the  coupled  linear 
system  (1.2),  (1.4),  (1.6)  and  (1.7).  Other  transmission  condition  (1.6)'  and  boundary 
condition  (1.7)'  can  be  treated  in  a  similar  way. 

Previously,  several  researchers  [1],  (2),  t3]»  [8]  have  studied  this  type  of 
problem  using  energy  identities  or  the  method  of  characteristics.  We  have  tried  both 
methods  for  system  (1.2),  (1.4),  (1.6)  and  (1.7),  yet  we  could  not  succeed  in 

establishing  an  affirmative  answer  to  [Q]  for  our  problem  for  any  wave  speeds  c,. 
c«  and  any  gain  constants  k*  ^0.  k*  a  0.  The  difficulty  probablv  can  be 

~0  “I  " 

inlerpretad  aa  followa: 

i)  At  the  coupling  point  x  :  1,  there  are  reflectiona  and  tranamieaiona  of  wavea 
which  the  energy  multipliers  are  not  sophisticated  enough  to  handle. 

ii)  For  coupled  linear  strings  the  method  of  charectorietics  worka  beat  when  the 
wave  speede  are  identical,  i.e.,  c,  :  c^,  cf.  [2],  yielding  sharp  decay  estimalee  for 
a  (cf.  (2.12)).  Otherwiaa,  this  method  is  not  convenient  for  coupled  strings. 

Of  course,  one  could  also  apply  the  nonlinear  Theorem  2  to  coupled  linear 
strings.  But  such  a  result  would  not  be  sharp  as  the  guin  coefricieiila  k'  and  k’ 

O  I 

becoae  severly  restricted  due  to  the  conditions  oh  the  spectral  radix  of  Bq  and  R,. 

A  recent  theorem  of  K.L.  Huang  offors  an  exlromly  UHofu)  way  to  prove 
exponentiiil  alobtlity  for  linoiir  aystema: 


Theorem  4  [7],  l.ot  A  bo  the  infinitoaimal  gonoralor  of  n  -annilgroup  oxpdA), 
which  satisfies 

lexp(tA)l  ^  Bo,  t  *  0.  for  some  Bq  ^  0.  (4.1) 


Then  exp(tA)  is  exponentially  stable  (i.e,  lexp(tA)l  <  Ke  K,a  ^  0,  t  *  0)  if  and 


c  fJ>  c  p(A),  the  resolvent  set  of  A; 


(4.2) 


emd 


Bt  ■  sup{|(iu  -  A)  *1  lu  c  R)  ^  • 


(4.3) 


are  satisfied. 


□ 


This  theorem  has  recently  been  applied  in  C4],  to  establish  an  exponential 

stability  result  for  an  Euler-BernoulH  beam  with  bending  moment  proportional 
control. 

Now  we  are  in  a  position  to  apply  Huang's  Theorem  4  to  (1.2),  (1.4),  (1.6),  (1.7). 
l.et  US  recast  the  problem  into  an  equivalent  hyperbolic  system:  For  x  e  (0,1), 
define 


u,(x,t)  •  5[-c,  l^(x.t)  +  |^(x,t)], 
Uj(x.t)  ■  ^[cj  |y<(2-x,t)  +  Jj^(2-x,t)]. 
u,(x.t)  •  |[-c,  |^(2-x.t)  +  |^(2-x,t)], 
U.(x,t)  ■  |[c, 

Than  (1.2)  becomaa 


U| 

-C, 

0 

0 

0 

u, 

i 

Ua 

0 

-c. 

0 

0 

i 

u. 

it 

Ua 

0 

0 

C» 

0 

ix 

u» 

0 

0 

0 

C,. 

u. 

Afler  Biraightforward  calculations,  we  gel  from  (1.4),  (1.6),  (1.7)  the  following 

boundary  conditions: 


At  X  =  0, 

u'(0,t) 

At  X  -  1 , 


u‘>{l,t)  =  4" 


c  ♦!(' 
I  o 


u>*(0.t)  «  DoU>'(0,l). 


( 1 .  r, ) 


2c,  -c,+c,+k,c,c, 

c,-c,*l<Jc,c,  2c, 


•  u*(l,t)  .  D,u*(l,t),  (4.6) 


where  4  =  Ci+Cj+kJciC,,  u'  =  (u,,u,}  and  u**  =  (u,,u»). 

Obviously,  the  underlying  Hilbert  space  is  X  •  (L'(0,1)]*  and 
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,  r't'i 


JH 


Pr 
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ti*»-  ♦ 


h  ;i  ^ 


|c  •♦^c  ■♦^k'c  c 

It  a  t  I  ai 


c  •♦c  ♦k^c  c 

1  a  I  I  a 


_  f(cj+Cidil£x£il’  +  lc.-cj*k?c.c.l* 

l|c  >c  ♦k^c  cl  Ic  ♦c  +k*c  c 
li  t  a  I  I  a/  (  «  a  I  I  aJ 


_  2  (c^*c,-)(,c,c,)(c,-c,*k,c,c.)  in 
(C.*c,.kjc.c.)* 


Bltfc.c* 


>k»c  c')«  ■  '=“•  f;)  ‘  °*  "  * 


(c  +0  +k'c  c 

I  >  I  I  a 


So  (4.14)  ia  varlfiod. 

Inequality  (4.15)  ia  obvious  from  the  definition  of  z,  in  (4.13),  because 


Ci-Cj-kiCiC, 

A 


c»“kiCir» 


Ic  *c  ♦k’c  c 
I  I  a  1  (  a  I 


<1,  for  any  k|  s  0. 


Therefore  the  proof  is  complete. 


Lemma  6.  The  operator  A  defined  in  (4.4),  (4.7)  satisrios  conditions  (4.1)  - 
(4.3),  provided  that  k'aO,  k*  fcO. 


Proof.  (4.1)  is  trivially  satisfied  because  A  is  a  dissipative  operator.  (4.2)  follows 
from  Proposition  5  because  equations  (4.12),  (4.10)  and  (4.9)  are  all  solvable  for  any 
given  f  e  X. 

Finally,  (4.3)  is  satisfied  because  in  (4.11) 


>VI 


IF.I  4  Klf  .1., 


j  =  1.2. 3, 4 


are  satisfied  and  because  of  Proposition  5,  (4.12).  (4.16).  (4.10)  and  (4.9). 


Hence  we  conclude 


Theorea  7.  Under  the  assumptions  k*  ^  0,  k*  ^  0,  the  energy  of  the  coupled  linear 
vibrating  strings  (1.2),  (1.4),  (1.6)  and  (1.7)  decays  uniformly  exponentially. 


Therefore,  for  coupled  linear  strings,  one  stabilizer  (kj  ^  0,  kj  =  0)  is 

sufficient  to  cause  exponential  decay  of  energy.  The  feedback  gain  k*  can  be 

0 

chosen  arbitrarily. 

By  straightforward  calculations,  one  can  easily  sea  from  (4.12)  that  Proposition 
5  in  general  does  not  hold  if  k*  =  0,  k*  ^  0  (i.e.,  only  one  stabilizer  is 

Installed  In  the  middle  of  the  span),  unless  C|  and  Cj  satisfy  certain  special 
relations,  such  as  C|/cj  =  2,  e.g.  More  importantly,  this  exponential  stabiliy  is  not 
robust  with  respect  to  Ci/Cj  in  the  sense  that  if  Ci/Cj  differs  just  slightly  from 
2  (or  certain  given  number),  then  exponential  stability  no  longer  holds. 

Therefore  we  see  that  a  point  stabilizer  installed  at  the  boundary  is  robust 
with  respect  to  exponential  stability  for  coupled  nonlinear  and  linear  vibrating 
strings. 
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